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ABSTRACT 
by 


Mario A. Baltazar 


Rational design methods based on a statistical 
description of both the ship load and strength have been 
proposed by several authors in lieu of the traditional 
bending moment calculations. In general, the effect of 
horizontal moments on the primary strength has been 
neglected. It is shown in this paper, as several others 
have found, that the horizontal moments can be very 
Significant. Two methods are described to obtain the 
magnitude of the equivalent combined moment. 


An example application is shown for a large tanker 
by use of the MIT seakeeping program. Wave moment trends 
with speed, sea state and wave heading are noted. Measured 
and predicted values of the vertical bending stress are 
compared. Finally, the parameters of the weibull distrib- 
ution to describe long-term wave moment amplitudes are 
derived. 
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Chapter I 
INTRODUCTION 


The traditional method for determining the load for 
primary strength calculations consists of poising the ship 
on a trochoidal wave of length equal to that of the ship and 
wave height equal to one-twentieth (1/20) of the length. 
The static bending moment is then determined by the load- 
shear-moment relations from beam theory. There have been 
various proposals to modify the wave height (e.g., H = Gael. 
or H = 1.1VL) partly to correct for dynamic effects and partly 
because it is thought than an L/20 wave is too severe for the 
very large ships. Notwithstanding, this empirically deiived 
procedure was considered entirely satisfactory until relatively 


recently and is still in prevalent use inasmuch as it gives 


conservative estimates. 


The appearance of large ocean-going tankers and bulk 
carriers and novel ship types gave impetus towards the search 
for a “rational" procedure for ship structural design. It 
was realized that the application of the standard bending 
moment calculations for these ships lead to unrealistic re- 
sults. At the same time new knowledge relating to ship 
structural design had accumulated and experimental and comp- 
utational facilities to advance this knowledge have improved 


tremendously. 





Many marine-oriented societies and organizations 
undertook basic research to improve the accuracy of pre- 
dictions of sea loads and motions, and the development of a 
rational design method applicable for all ship types. Full 
scale and model tests data were collected. The random 
nature of both the loads on a ship and its capacity to re- 
sist these loads was recognized, so that probability theory 
was utilized to develop statistical models for the design of 
ships [1-14]*. The advantages of a probabilistic approach 
over the traditional deterministic calculations are discussed 


imsreterence [7]. 


In the conventional method, the primary strength of 
a ship is based solely on the longitudinal stresses caused 
by vertical bending, or bending about the horizontal neutral 
axis. The effects of lateral bending, or bending about the 
vertical neutral axıs, is neglected. However, for ships 
travelling oblique to the waves, the lateral loads reach a 
Magnitude equal to or greater than the vertical bending 
moment. The combined effect of lateral and longitudinal 


bending may be critical for some ship types. 


A requisite for an integrated probabilistic design 


procedure is a more precise knowledge of the loads acting on 


*The numbers in brackets refer co references listed 
in the Appendix 





a ship. This paper represents a small portion of the total 
effort to completely define the structural loads experienced 

by a ship in a seaway. It is an investigation of the magnitude 
of combined horizontal and vertical bending in irreguiar seas. 
Due to their relative phase difference, the amplitudes of the 
vertical and horizontal moments generally are not acting at 


the same time. 


In the analysis, the loads are treated as random 
variables. Consideration of the total longitudinal stress at 
the deck leads to the definition of an equivalent combined 
bending moment. Two methods are developed to generate the 
Statistical parameters of the combined load so that the usual 
statistical methods to predict long-term extreme values of the 
load can be applied. An example application is then made for 
a large tanker which had been previously instrumented so that 
full scale data for the vertical bending stresses experienced 
in service are available for purposes of comparison. Such a 
comparison serves to check yet another time the validity of 
computational techniques to simulate ship motion and loads in 
waves, and in particular the MIT 5-D seakeeping program. The 


predicted moment responses are compared and analyzed for trends. 





Chapter II 


PROCEDURE 

II.l Determination of the Equivalent Combined Wave Bending 
Moment 

A ship travelling at an angle to the waves will be 
subjected to unsymmetrical bending such that the stress at 
the deck is not uniform but assumes a maximum value at one 
edge, as shown in figure II.l. This primary stress may be 
thought of as the combined effect of vertical and horizontal 
bending and, to a much lesser extent, torsion. Previous 
studies [15, 16] have established that torsional moments are 
relatively small, never exceeding 10$ of the vertical mcment 
in any given condition. In this study the effects due to 
torsion are considered negligible, but it is recognized that 
torsion loading may be critical for some ships, as those 
with large hatch openings [17]. The contribution to the 
ship primary stress may be included, if so desired, in a 


manner analogous to the method described here [18]. 


The stress at the deck edge is given by 








Vertical Bending Stress Horizontal Bending Stress 
(sagging) 





Combined Stress 


Figure II.1 


The Longitudinal Stress Distribution 
on a Ship Hull Girder 





where 
On = combined, or total, stress 
Sy = stress due to vertical bending 
I,” stress due to horizontal bending 


In terms of moments, one can write 


Q 

I 
AF 

+ 
au 


where 


My = vertical bending moment 


M = horizontal bending moment 
_ Iyy . . 
ky . a vertical section modulus 
Z 2 : : | 
ku = : xd horizontal section modulus. 


multiplying both sides of the equation by Ky one gets 


ky 
Onky = M, + ko M 





If one defines an equivalent combined bending moment 
M = o Kv and the section modulus ratio K = —, then one 


y 


gets the relation 


Mp = My + My (1) 


It is emphasized here that one cannot simply add 
the maximum values of the vertical and horizontal moments to 
get the maximum combined effect. In general the moments at 
any point along the length of the ship are not acting in 
phase. Furthermore, it may be of interest to specify the 
complete distribution of the primary stress at the deck. 
Within the limits of the beam theory, this is done by getting 
the stress at one other point across the ship's deck. A 
convenient reference point is at the ship's centerline where 
the stress due to lateral bending is zero for most ships due 
to port-and-starboard symmetry, i.e., Om = Oy at the centerline. 

In the statistical approach to sea loads, the ex- 
citing force is assumed to be a stationary, narrow-band and 
gaussian process with zero mean over a short period of time 
(e.g. 30 minutes). The ship is considered to be a linear 
system so that the response is also stationary, narrow-band 
and gaussian with zero mean over the same period of time. 
Then the amplitudes follow the Rayleigh distribution [19, 20]. 
Over a long period of time, the assumption of a stationary 


excitation process cannot be made. However, it was found 





that the amplitudes follow closely the Weibull distribution 
[8, 9, 10J. It is shown that the Rayleigh is a special case 


of the Weibull distribution given by: 
E s-1 8 
f (x) = (SKIK) exp -(x/k) X20 (2) 
FY, (x) = 1 - exp -(x/k) > X20 (3) 


wher f, (x) and F. (x) are the probability density and distrib- 
ution functions respectively, and s and k are parameters. To 
get the Rayleigh distribution, one sets s=2 and k = /E where 
E is the mean square value of the amplitude of the response. 
Furthermore it was found that the best-fit Weibull curve is 
approximately the exponential distribution obtained from 

Aqua tion (2) and (3) by setting s = l and k = m where m is 
the mean value of the response amplitude. The statistics of 
the ship bending moment response, therefore, can be specified 
if the mean square value, or variance, of the process is 


known. 


The usual method to obtain the variance of the re- 
sponse in random seas, attributed to St. Denis and Pierson, 
consists of specifying the energy density spectra of the 
Ocean waves and obtaining the ship respcnse amplitude opera- 
tor (RAO), or the response per unit amplitude of regular 
wave, by model tests or some computational procedure using 
the strip theory approach. From the theory of linear systems, 


the spectral density of the response cen be found from the 





relation [21]: 


Sg — | H(w) | 7S (w) 


where Sp (4) = spectral density of the response 


H (w) 


ship response amplitude operator (RAO) 


S (w) 


Spectral density of the seaway. 


The area under the response spectrum is the variance 
of the process and represents the sum of the responses to the 
infinitely many sine waves that are thought to compose the 


actual sea. 


Two methods used to obtain the variance of the wave 
bending moments by computer simulation techniques are ex- 


plained below. 
1) First Method 


Equation (1) can be written showing the moments as 


a continuous random function of time, t: 


Mn(t) = M,(t) + KM, (t) 





10 


Using the assumption that the sea can be represented by in- 


finitely many harmonic waves of small amplitudes, continuous 


frequency and random phase, one can write the vertical and 


horizontal moment responses as: 


where 


then 


Mt) 


Or 


Mp (t) 


which can be 


M; (t) 


V cos(0 t + € ), and 
n n n 


ME) = H Cos (o t + 0) 


amplitude of the vertical moment response to the 


nth wave 


amplitude of the horizontal moment response to 


nth wave 
frequency of the nth wave 


phase angle of the vertical moment with respect 


to a wave with its crest amidship 


phase angle of the horizontal moment with respect 


to a wave with its crest amidship 


oo co 


er V cos (wt + E) + K 2H cos(u t + 9 ) 


on 


* e | 
„L [v cos(u e + e) + KH cos (v t + i 


rewritten as 





PL 


oO 
= + 
M(t) 2 [(v,cos e, KH, cos 0. ) cos wit 


- (V sin €. * KH sin 0 ) sin o d 
n n n n n 


Let 
X = V cos € + KH cos Ü 
n n n n n 
Y = V sin € + KH sin Ü 
n n n n n 
Then, 
co 
MC) — [x cos w t + Y sin o d (4) 
e n n n n 


Since Mq is also a continuous random function of time, it 


can also be represented by: 


l| 0-38 


Matti = 
n 


I ne + 0. ) (5) 


A comparison of equations (4) and (5) shows that 


-1 Y 
0 = tan n, and K = X” + y? 
n = n n n 





12 
Then 


2 = V 2 + vin $ + 2KH V cos(e - 0 ) (6) 
n nn n n 

Equation (6) gives the amplitude of the combined 
moment due to vertical and horizontal bending, with the 
phase difference in the occurence of the component peaks 
taken into consideration. The MIT 5-D seakeeping program 
calculates the ship motions and the vertical and horizontal 
moment responses to regular waves (RAO). A slight modifica- 
tion of the program to obtain the equivalent combined moment 
response operator was made by use of equation (6). The RAO's 
were then used in the usual power spectral analysis in the 
frequency domain to obtain the variance of the responses in 
irregular seas. Comparison of predicted and experimental re- 
sults have shown generally good agreement [15, 16, 21, 22, 23] 
and have led to widespread acceptance of the superposition 


method of St. Denis and Pierson. 


2) Second Method 
Considering the wave loads as random variables, it is 


generally accepted that the bending moments My and M are sta- 


tionary, narrow-band, gaussian processes with zero means for 


short periods of time. Then Mn is also stationary, narrow-band 


and gausSian with zero mean, i.e., in statistical notations: 


E Mn = EM, + KE My = 0 





Jes 


Furthermore, 
E m3] = = |( + 101,9” 
p[m2| = [ve] ` ele | + 2KE [My - M 
but the covariance is defined as: 
s vens 5 Po? | (My n a, Viel, E ds 
where My = m. = mean values = 0 


ler = correlation coefficient of the vertical and 


horizontal bending moments 


therefore, 


déi = spe] + Sell + 
2% e e |n |) (7) 


For a zero-mean process, the mean square is equal to the 


variance. 
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As in the first method, the value of the section 
modulus ratio, K, must be known for a given ship or esti- 
mated during preliminary design. An attempt was made to 
examine the section modulus ratio of several sizes of tankers 
to determine if there was a trend with the growth in the size 
of these ships. As shown in figure II.2, the points are 
widely scattered with no apparent trend. Since the value of 
K depends on the construction of the ship (framing system, 
Scantlings, etc), perhaps the scatter of points is indica- 
tive of the present uncertainty about both sea loads and the 
strength of ships, and the rather arbitrary way a "safety 


factor" iS arrived at. 


euse, Che value of the correlation coefficient, 
fug! of the vertical and horizontal moments must be deter- 
mined.  Implicit in the use of equation (7) is the assumption 
that an average value of Pug already takes into account the 
phase difference of the peaks of M, and M in all conditions 


at sea. 


With this method, the bending moment response ampli- 
tude operators may be determined by experiment or computer 
techniques. Spectral analysis could then be performed to 


obtain the variances of My and M and finally the combinea 


H’ 


moment Mo. 





ection Modulus 


DWT(X10? Tons) 
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Figure II.2 
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Section Modulus Ratio VS DWT of Tankers 
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II.2 Long-Term Wave Load Analysis 

After the statistics of the short-term response are 
determined for the range of sea states that make up the 
actual seaway, the expected long-term loads may be predicted. 
Several methods to obtain the long-term distribution of wave 
loads and the prediction of extreme values have been proposed 
(e, 9, 11, 18, 231]. Mansour (18] has described procedures 
for the analysis of failure under extreme loads under two 


conditions: 


1. Short-term analysis - when the route of the ship 
is not known and a "severest condition" criterion 


is imposed. 


2. Long-term analysis - If the ship's route is known 


and if that route is more or less permanent. 


In this study, a procedure is discussed to determine 
the value of the parameter k = m of equations (2) and (3) 
for a ship travelling oblique to waves in all sea conditions 
described by the significant wave height and period (h,T) as 


in [24]. The following assumptions are made: 


l. wave headings are uniformly distributed over all 


directions. 





17 
2. The speed is essentially determined by the state 
of the sea. The speed-significant wave height relation is 
adopted from [14] as follows: 


Speed, V knots significant wave height, h ft. 





16 less than 12.5 
12 posto 21.5 
6 greater than 21.5 


The procedure followed is described below: 


l. Obtain the RMS response to all sea states (i.e., 
range of sea spectra) in all conditions of speed 
and wave headings by some analytical method as 


the "MIT 5-D seakeeping program." 


2. For each speed, get the average of: the RMS 


response over all wave directions. 


3. Use Bretschneider's two - parameter family of sea 
spectra to get the response to any sea described 
by a particular (h,T) pair. This is most easily 
done by scaling the responses to the sea spectra 
of step (l). We call this matrix of the response 


R with elements Ze 


4. For each wave height in [24], determine the 


probability (i.e., relative frequency of occurence) 
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of each observed period. We call the resulting 
matrix of probabilities P with elements ZEE 
Determine the average RMS response for each 


significant height by the relation 


where the subscripts i and j refer to rows and 
columns of a matrix. The resulting column 


matrix is called M. 


Determine the probability of occurrence of each 
significant wave height. We call this column 


matrix Q with elements d. 


The mean or expected value of the RMS of the 


response for a particular area is given by 


mean RMS = ` m. x ae 


Then the desired parameter m = 1.25 x (mean RMS 
of the response) in the case of a short-term 
analysis [18]. For the long-term analysis, a 


time factor is additionally introduced to get the 





ÉS 


weighted mean of the RMS of the response for all 


areas in the ship's route. 


Steps 4-6 is equivalent to getting the probability 
of occurrence of each (h,T) pair in any one area being de- 
scribed. If we call this matrix P' with elements Page then 

1 


RMS = T Me 
mean } e a x PO 


The data in [24] have been computerized and could 


be used in the above calculations. 
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II.3 Details of the Procedure 
The methods for wave load analysis and prediction, as 

explained in the foregoing sections of this chapter, were 
applied in the analysis of a large tanker. The UNIVERSE 
IRELAND was one of five vessels of different sizes instrumented 
for full-scale measurement of stresses experienced in actual 
service. However, the strain gage arrangement on deck was 
such that the stress due to lateral bending was cancelled out, 
so that only vertical bending stresses were reccrded. There- 
fore, the results derived for the response of primary interest, 
the combined moment, could not be verified. The description 
and analysis of the results of the 3%-year project sponsored 

by the American Bureau of Shipping (ABS) is reported in [25] 
and provide the basis for comparison with the analytically 


derived values of vertical kending stress at amidship. 


The general characteristics of the UNIVERSE IRELAND 


are listed below: 


all aras placement, tons... +... < < < | 325548117 
rana ab overall, ft a e 1135.17 
Length between perpendiculars, ft . . . 1083 
buco ge ed is 174.87 
Demian o 105. 
Design draft (keel), ft . . . +. +< < < ° 812025 
BeelBesetfricient (wl). . + » 2... 0.86 
Section modulus, top, in?2-ft . . . . . . 566,794. 
Section modulus, side, in2-ft . . . .. 617,350. 


The ship's lines and other data pertaining to the 


ship nece.sary to gcnerate the input for the computer analysis 





21 


were obtained from the ABS. It was necessary to use an 
approximate lightship weight distribution because of an 
apparent discrepancy in the distribution curve provided. 
Since the still water bending moment, at least, is very 
sensitive to the distribution of the weight, a probable 
error was anticipated in the calculated values. For ex- 
ample, the maximum still water bending moment calculated 
by the computer exceeded the value determined by the ABS 


by about 11.5%. 


The MIT 5-D seakeeping program is fully described 
in [26]. Only minor modifications were necessary to do the 
quasi-static combined bending moment analysis. The re- 
sponse to long-crested irregular seas were used in conjunc- 
tion with an auxiliary program to obtain the responses to 
short-crested seas. This program uses the cosine-squared 


spreading function 
2 2 
SA) = 7 cos (+) 
where “= angle between the wind direction and the wave 
components. The principal wind directions were given ex- 


actly the same values as the wave-to-course headings. 


The computer was run for ship's speed of 6, 12 and 


16 knots, nine wave headings from 0° to 180° in increments of 
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22,5, and nine sea states described by the significant 


height and period. The computer output consisted of the 
root mean square of the response, the significant height 

(u É) and the average of the 1/10th highest (g A?) of the 
response measured peak-to-peak for all combinations of the 
three speeds, nine wave directions and nine sea states. 
Since it was desired to compare the values of the vertical, 
horizontal and combined wave moments, output was obtained 
for all three responses. Furthermore, to be as accurate 

as possible, the three EON ss were calculated for twenty 
stations specified along the length of the ship, even if 
eventually only the responses amidships were utilized in 


the analysis. 


The sea data were carefully chosen to match the 
actual seaway experienced by the ship during the instru- 
mentation period. This was done by obtaining sea data from 
[24] for the areas lying along the service route of the 
UNIVERSE IRELAND, which made regular runs between the Persian 
Gulf and Western Europe by way of the cape of Good Hope [25]. 
Although the seakeeping program can accept and use an actual 
sea spectrum or simulate developing or decaying seas by the 
Bretschneider two-parameter spectrum, the Pierson-Moskowitz 
one-parameter family for fully developed seas was chosen. 
Nine sea states with significant heights ranging from 3.6 feet 


to 43.3 feet (corresponding to observed periods of 1.5 seconds 
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to 18.5 seconds) were used. Spectral ordinates for forty 
frequencies were calculated while twenty five points were 
used to describe the response amplitude operators covering 
the wave frequency range 0.1 < SÉ 9.0 with most of the 
points spaced closely around the expected spectral peak at 
SE 1.0. A is the wavelength and L is the length of the 


ship. The response to regular waves (RAO's) were not 


printed in order to save on the cost of running the computer. 





Chapter III 
RESULTS 
The following general areas of study were explored: 


l. Comparison of the derived results with those 


in other published works on the subject. 


2. Comparison of the variance of the combined 
moment obtained from the two methods described 


in the previous chapter. 
3. Comparison of theoretical calculations for the 
vertical bending stress with full-scale 


measurements. 


4. Examination of bending moments trends and wave 


load prediction. 


The results of the above studies are presented in 


the following sections. 
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III.1 Comparison of Results with Other Published Works 

The calculated value of the root-mean-square (RMS) 
of the bending moment responses are plotted in figures III.1 
to III.4 as a function of the wave headings. For contrast 
this was done for the responses to long-crested and short- 
crested fully developed seas with significant wave heights 
Dr 7.4 ft., 19.2 ft. and 43.3 ft. It was intended to show 
the response to a low, moderate and high sea states. Since 
only the value of the horizontal moment in "realistic" seas 
are of interest, the RMS of the response were plotted only 
for short-crested seas. The graphs were prepared for a 
representative ship's speed. The general behavior of the 


response did not change with speed. 


Faltinsen [22] describes the comparison between 
theory and experiments of tne wave-induced horizontal and 
vertical forces and moments of a series 60 hull form with a 
block coefficient of 0.80 and a Froude member of 0.15 moving 
in regular waves with different headings at an angle to the 
direction of the waves. The predictions made by the new 
strip theory developed by Salvesen, Tuck, and Faltinsen, 
which is used in the MIT 5-D seakeeping program proved very 
satisfactory. The theory is described by Salvesen et al, in 
[27]. For vertical bending moments, the theory predicted 
higher values than the experiments for ° < 0.7 in quartering 


and following waves and lower values for A/L < 0.6 in bow 


waves. This trend may also be observed by comparing the RMS 
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values of My in short-crested seas for wave headings of o» 
(following seas) and 180 (head seas) in figures III.2, 3 
and 4. For horizontal bending moments, appreciable diff- 
erence existed between theory and experiments only for beam 


Seas. 


One of the early investigations of the effect of 
horizontal moments in oblique seas is described by Numata 
[16]. It was established in this study that lateral moments 
reach values that are substantial and that peak bending 
moment exists in the region of the speed for synchronous 
pitching and heaving motions.  Wahab [15] shows that the 
importance of the horizontal bending moment increases with 
increasing ship length. For the sea condition used in his 
investigation, the significant midship horizontal moment for 
a ship of 200 meter-length was 46 per cent of the significant 
vertical moment while for a 400 meter-length ship the per- 
centage increased to 69. The concept of effective wavelength, 
Los described by Numata was helpful to explain the results 
obtained from the computer analysis. However, only one sea 
spectrum was used in his study and some conclusions were 
reached that are not necessarily true for all conditions at 
sea for any given ship. Inasmuch as the RAO of the UNIVERSE 
IRELAND was not recorded by the computer, the figures in 
[15] proved very useful. The ship studied in that work is 
the same as the model described by Faltinsen and has general 


characteristics similar to the UNIVERSE IRELAND. The respense 








of the model to regular waves could be assumed 
stantially identical in shape and trends. The 
of amidship forces and bending moments in [15] 


use of a single spectrum representing a severe 


3l 
to be sub- 
calculation 
also made 


sea. 
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III.2 Comparison of the Calculated Values of the Combined 
Moment Response 


Table III-l is a summary of the results obtained 
for the RMS value of the combined moment response for the 
vessels' service speeds of 16 and 6 knots. The first set 
of values are those calculated from the RAO developed by the 
computer by use of equation (6) of the first method. The 
second set of values were calculated from equation (7) of 
the second method using the values en = 0.32 obtained from 
[18] as reprinted from the published result of a Det Norske 
Veritas Research Report. The reprinted source is shown as 
figure III.5. The values of fon represent the correlation 
coefficient of the most probable largest values of the 
vertical and horizontal moment for one ship life of about 
twenty years in the North Atlantic. They were meant to be 
reference values and not for use in design. The difference 
in M, values is shown for each sea state. It is immediately 
apparent that the results obtained by both methods are not 
very different. The third set of values was calculated also 
by use of equation (7) but using a value E obtained by 
averaging the responses in short-crested seas by the first 
method for all wind-wave directions and for TOS represent- 
ative spectra for low, moderate and high seas. The mean 
value on „tained in this manner was 0.532, significantly 
higher than the ISSC value. The difference in the results for 


Mn is also indicated for comparison. 
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Series 60, All Headings Included, Q=10 Š 
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Figure III.5 


Correlation Coefficient Between Vertical Bending Moment , 
and Horizontal Bending Moment 
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Lit. 3 Comparison of the Calculated Vertical Bending Stress 


and Full-scale measurements: 

Strain gages were installed on the port and starboard 
sides at amidship and several other stations along the length 
of the UNIVERSE IRELAND to meaSure stresses on the deck. The 
arrangement was such that only vertical bending stresses were 
measured. Therefore, though the combined moment could not be 
verified, a comparison between full scale data and the com- 
puter-derived values was possible. This comparison by 
weather group as specified by wind measurements (Beaufort 
scale) is shown in figure III.6. A summary is also pre- 
sented in table III.2. The conversion from bending moment 
to stress is done by the relation: Bending Moment = Stress x 
Section Modulus. The values compared are the average RMS 
values of the peak-to-peak stress for each wind scale. The 
full scale data were read off figure 11 of [25]. The con- 
version from wind to wave was done by the graphical relation 
shown in figure III.7 [23]. This conversion is summarized 


Der table III.3. 


The use of nine sea spectra with each significant 
wave height enables the calculation of the standard devia- 
tion of the predicted values of the vertical stress. For 
comparison with the measured data, these standard deviations 
had to be modified to obtain values based on wind measure- 


ments. This was accomplished by the relation [23]: 
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2 
S um S — + 1S a AH [tan 9 
1 2 3 == 
2 
where S, = standard deviation of the response with respect 
to the wind 
S, = standard deviation of the response with respect 
to waves 
S, = standard deviation of the wave height with re- 


spect to the wind, as 


from 1231. 


AH = wave height increment 


shown in figure III.8 


6 = average Slope of the response vs. significant 


height plotted in figure III.9. 


For comparison, table III.4 of the 


all the cases studied is drawn up. 


Standard deviations for 
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Relationship Between Significant Wave Height 
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III.4 Examination of Bending Moment Trends and Wave Load 


Prediction 

Since the traditional primary strength calculations use 
the vertical bending moment in a long-crested sea of arbirtary 
height as a basis for ship load, it is of interest to compare 
the values of the vertical bending response in long-crested head 
seas with the combined moment in short-crested seas for the 
higher range of the wave significant heights. Figures III.10 
and III.11 are drawn for sea spectra representing fully developed 
moderate and severe seas, respectively moderate seas, the com- 
bined moment for head seas exceeded the long-crested vertical 
moment by about thirty sever per cent (37%) and in severe seas, 


the difference was sixteen per cent (16%). 


Similarly, it was deemed desirable to compare the rel- 
ative values of the vertical and horizontal moments for all 
conditions at sea. This is done in figures III 2, 4 and 6 for 
the responses in short-crested seas. One can see that the 
horizontal moment is not small compared to vertical moment; 
and in low to moderate seas, the former is at least as great as 
the latter. The above mentioned comparison of the response in 
long-crested and short-crested head seas reflect the relative 
contributions of the vertical and horizontal moments to the 


total load. 


Table III.5 illustrates the change with speed in the 
bending moment response for low, moderate and high seas. The 
tabulated values are the mean RMS peak-to-mean responses in 


short-crested seas averaged over all wind-wave directions. 
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Comparison of Vertical Moment Response in Long-Crested Seas VS the 


19 2 EE 


Combined Moment in Short-Crested Seas for h 
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In order to illustrate the procedure to obtain the 


parameter k for the long-term distribution of the wave re- 


sponse, a "Short-term" analysis was done for the worst 


condition in the ship's route (area 7 in [24]). This is the 


North Atlantic in the months of December to February and the 


tabulated data of [24] were used. 


given as tables III.6 to III 


this analysis, the result is 


mean RMS 


Since this is a "short-term" 


m = 1.25 (mean RMS) = 459775 


The results obtained are 


.9 for the combined moment. For 


367820 ft-tons 


analysis, one can then say 


fE. CONDS- 


The equivalent values for the vertical and horizon- 


tal moments were also computed. 


mean RMS response, My 


mean RMS response, M 


H 


These are: 


I 


197307 ft-tons 


231043 ft-tons 


Note that the expected value for the horizontal moment is 


larger because the ship encounters more low-to-moderate 


seas in any iength of time than severe seas. 








The mean value of the RMS of the response was 
calculated for all the areas in the route of the ship, 


a "long-term" analysis was conducted. Both My and M 


calculated and the results are summarized in table IIi. 


also 
i.e., 
were 


TOS 





Ki 
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Table III.6 


Mean RMS of the Combined Moment Response in Short-Crested Seas 
Averaged Over Uniformly Distributed Wind-Wave Directions 
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Table 111.10 


Mean RMS Value of the Response for all Areas in the 


Route of the UNIVERSE IRELAND 


Area | ELRMS Jof M, Response ELRMS Jof M, Response Time Factor 


154254 
120713 
109398 

66961 
110271 

80454 

92992 
108754 
146576 
123456 
161230 


Time Weighted Average: 


(RMS), = 124,562 ft-tons 


(RMS) = 517,450 ft-tons 


299294 
246724 
218546 
138544 
220161 
170798 
189706 
203516 
282585 
246480 
260578 


Ñ N N FP FP a a | PF IM t 


= 155,/02 ft-tons 


ES l 


m, = 646,813 ft-tons 


= 4.15 


Ef! 








Chapter IV 


DISCUSSION OF RESULTS 


IV.1 Evaluation of the analytical Results 


The results obtained in this study agree with the 
previous findings in other studies that the lateral moment 
can be very significant for ships in waves. Some definite 
trends of the moment responses can be observed from figures 
III.1 to III.4. From the results in short-crested seas 
shown in figures III.2, 3 and 4, one sees that horizontal 
and vertical moments change very little with wave direction 
in moderate seas. On the other hand, there is a Dbivad peak 
away from the head and following seas for all moment responses 
in low sea states. In severe seas, the maximum response for 


My and Mn are in the head and following seas while M. stays 


H 
almost constant regardless of wave heading. For the long 
crested seas of figure III.l, one can see that the vertical 
moment peaks at about 60% for low to moderate sea states and 
assumes extreme values at head and following seas in high 
seas. The equivalent combined moment has two peaks in all 
cases. The peaks occur in wave headings near 60° and 120° 


with a shift away from the beam seas case as the sea grows 


more Severe. 
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One observes from figure III.2 and III.3 that the 
horizontal moments exceed the vertical moments in almost all 
headings for low and moderate seas. However, in severe seas, 
for wnich structural design is based, one sees the customary 
relation of the vertical moment exceeding the horizontal 
moment, especially in head and following seas. The reason 
for this phenomenon is that lateral bending is more respon- 
Sive to higher frequency excitation than vertical bending. 
One can see from figures IV.l and IV.2 that the spectral peak 
of the lateral moment response for all wave headings are gen- 
erally situated at higher frequencies than the corresponding 
response curves for the vertical moments. One can note fur- 
ther that the spectral ordinates of the RAO for wave headings 
00 , 110° and Eer are higher for the lateral moment than the 
vertical moment. Therefore, when the energy of the exciting 
waves are concentrated in the higher frequency range (i.e., 
in low sea states), the lateral moment exceed the vertical 


moment. In high sea states, the reverse trend is true. 


For a similar reason the peaks of the response curves 
for the vertical moment occur away from head and following 
seas in low sea states. The shift in the location of the 
spectral peak of the RAO's in the direction of higher cir- 
cular frequency, w, can be explained by considering the 
effective wavelength in oblique seas. It is generally known 


that the wavelength nearly equal to the ship length causes 
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Vuave 


Figure 1V.3 


The Effective Wavelength in Oblique Waves 
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the maximum midship bending moments. In oblique seas the 
ship senses the effective wavelength, La’ as shown in 


figure IV.3. The relation with the actual wavelength is 


Lo = L/cos ß 
where Lo = effective wavelength 
L = actual wavelength 


^0» 
Il 


heading angle 


AS P approaches 90? from either 0? or 180 the actual 
wavelength decreases while the value of Lo stays fixed. 
Therefore the peak frequency of the response goes up, thereby 
shifting the RAO spectrum nearer the sea spectral energy 


density peak in low sea states. 


The theoretical calculations predict greater vertical 
moment response in following seas than in head seas, contrary 
to experimental results. Wahab's results [15], shown as 
figure IV.4, confirm this apparent discrepancy between the 


theory and experiments. 


Finally, one notes that the combined moment in 


"realistic" seas is significantly higher than either of the 





components in all cases but, because of the relative phase 


angles, is not equal to the sum of their individual peaks. 
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IV.2 Comparison of the Values of the Combined Moment 


Table III.1 shows the comparison of the response 


RMS value of the combined moment, M as obtained by the two 


mi 
methods described in chapter two. It can be readily seen 
that the results are not Significantly different, the aver- 
age for all sea spectra being on the order of eight per cent 
(8%). This confirms the statement that the average value of 
the correlation coefficient, Ar already considers the phase 
difference of the vertical and horizontal bending moment 
peaks. As would be expected, the third set of values for M, 
were closer to the computer values. It is realized that the 
two values c that were used were derived in a different 
manner and perhaps should not be directly compared. In de- 
riving the second value, it was found that the cociiicient 
varies strongly with the state of the sea and the wave di- 
rection and therefore that no single value will be "true" for 
all conditions at sea. If this method is to be adopted, 


there is need to study the correlation coefficient further. 


There 1S a dearth of published material on the subject. 


The calculated value of the combined moment also de- 
pends on the relative phase angle of the vertical and hori- 
zontal moments. Therefore, it 1S necessary that the pre- 
dicted values be at least as accurate as the predictions for 
the vertical and horizontal moment responses. Comparison of 


theoretical and experimental phase angles shows that the 
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agreement is in general quite good [22]. The maximum dis- 

crepancies occurred for wave headings of 50° and 110° £or 

the wave lengths A/L = 1.2 and ` = 1.5. (Figures IV.6 and IV.7) 
In view of the above mentioned considerations it 

would seem to be indicated that both method 1 and method 2 

to obtain the combined moment are equally "good". The first 

method is probably more convenient for use in analytical 


studies. 
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Bl 


67 





nn QO ----- Experiments, Wahab 68 
DEL go. lk. Theory 





Wave Directions 110° 


(Head Sea = 180 ) 
£- 9 


2526 231.0 





-100 


Figure IV.7 


Phase Angle Between Horizontal and Vertical 
Wave Bending Moment at Midship 
for Series 60, Block 0.80 
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IV.3 Comparison Study of the Vertical Bending Stress 


Data on the full-scale measurement of vertical 


bending stresses on the deck of the UNIVERSE IRELAND can be 


ound in [25]. 


A comparison between calculated and measured 


values was done by plotting the average RMS peak-to-peak 


amidship bending stresses. The analytical values were ob- 


tained by calculating the response to each of nine sea 


spectra corresponding to the observed periods in the tabu- 


lated data of [24]. The significant wave heights are the 


wind-to-wave conversion shown in table III.3. The average 


RMS peak-to-peak response is the weighted mean of the cal- 


culated response to each sea specified by the significant 


height. This is best described by the relation 


where 


NX 


DAT) 


E, > p(T) 


mean RMS V 


lI 
Il AO 


RMS peak-to-peak bending stress 


relative frequency of occurrence of each 
observed period for a given significant 


wave height. 
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The wind-wave relation as described by the Pierson- 
Moskowitz random and ITTC 1966 line, the corrected Roll curve 
for the North Atlantic and the corrected Yamanouchi curve for 
the North Pacific were examined. The conversion is summar- 


zed in table III.3. 


The use of the Pierson-Moskowitz Random line pro- 
duced the worst case in the comparison study. While the 
wind-to-wave relation is intended to convert wind measure- 
ments to significant wave heights for seas in random stages 
of partial development, there is a bias in favor of seas that 
are near full development [25]. Consequently, the values 
based on this wind-to-wave conversion, coupled with the known 
tendency of the theory to over-predict [21], are inuach higher 
than the measured stresses. The Roll and Yamanouchi curves 
were based on observed wave heights and corrected for signif- 
icant heights. They are believed to be better representation 
of an actual sea. Since the ship route is neither entirely 
in the North Atlantic nor in the North Pacific, one sees fairly 
good agreement with results based on both the Roll and 
Yamanouchi relations for the lower half of the sea states as 
described by the beaufort scale. It can be seen that theoret- 
ical results over-predicted the average values of the RMS 
stress, specially in the higher sea states. There are several 


probable causes of this: 
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Since the largest divergence of compared values are 
consistently in the more severe seas, one is led 
to believe that the assumption of linear motions 
for the ship, in particular resonance in roll, 
exaggerates the analytical values. Obviously, 
heave and pitch are large in heavy seas. The 
divergence caused by roll resonance gets bigger 
as the energy of the waves are concentrated more 
and more in the lower frequencies. Faltinsen [22] 
used a first approximation of roll at resonance 
aS input values. It is also to be noted that 
the ship was not provided with bilge keels in 
the computer analysis. It is not known if the 


acutal ship has these dampeners. 


Loukakis [21] suggested that deflections of the 
hull can relieve some bending stresses, so that 
even if bending moment predictions are accurate 
the resulting conversion into stresses may result 


in over-prediction. 


The inexact weight load distribution used in the 
analysis could have caused some of the over predic- 


tıon. A qualitative measure of this effect can be 





12 


gained by the excess in value in the calmer seas 
when the contribution of inertia effects to the 


dynamic loads are relatively greater. 


4.  Exaggeration of the reported occurrence of the 
higher period waves. In high seas, when the ship 
is pitching and heaving, it is understandable 
that such errors could be made by untrained ob- 
servers. In this study, this effect would be 


ME minor. 


5. Inaccuracy in the measured data due to insuffi- 


cient number of records in extreme weather. 


For direct comparison of the scatter of the derived 
values of the stress with the measured values, it was necess- 
ary to determine the standard deviations with respect to 
wind measurements. One can see from table III.4 that the 
calculated response to the nine sea spectra with the same 
significant wave height were more widely scattered than the 


measured RMS peak-to-peak values. 





73 


IV.4 Bending Moment Trends and Wave Load Prediction 


One can observe from figures III.10 and III.11 that 
the combined moment in short-crested seas are considerably 
higrer than the vertical moment in head seas. The difference 
in severe seaS 1S in the order of 15-20 per cent. Therefore, 
at least for this ship, calculations based on vertical bend- 


ing only could lead to disastrous results. 


Comparison of the horizontal and vertical bending 
moments in short-crested seas have been partially discussed 
in section IV.1l. It was explained that in low-to-moderate 
seas the lateral bending response to wave excitation is 
greater than vertical bending. One consequence of this is 
its possible significance in fatigue failures. nother is 
that over any length of time the expected wave bending load 
will be greater for lateral moments than for vertical moments. 
This 1S so because the ship will encounter more low-to-moderate 
seas than high seas in its lifetime. The result for the 
short-term analysis for the North Atlantic confirms this 
statement. The statistics reported in [24] were collected 


for a period of seven years. 


Table III.5 shows that in short-crested seas the 
bending moment varies only slightly with speed. This state- 
ment holds true for all the moment responses, with My 


having the greatest variation for all three sea states shown. 





This indicates that reducing speed in a seaway reduces 
dynamic stresses due toslamming, etc. but will not 
Significantly reduce the primary stresses due to wave 


bending loads. 
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Chapter V 


CONCLUSIONS AND RECOMMENDATIONS 


The subject of the combined effect of horizontal and 


vertical wave moments has been investigated. A number of 


conclusions can be drawn from the results of the study. 


Eier 


The results obtained by the MIT Seakeeping 
Program agreed well with published results of 
theoretical and experimental investigations 
of ship motions and loads in a seaway. The 


following bending moment trends were „userved: 


The horizontal moments experienced by the UNIVERSE 
IRELAND are not negligible in comparison with the 
vertical moments in almost all conditions at sea. 
In general, the relative significance of the lat- 
eral moments increase with the size of the ships 
[15]. In low sea states, the lateral moment 


exceed vertical moments. 


The variation with wave headings of the moment 
response seem to depend on the state of the sea. 
For long-crested seas, the maximum vertical 
moment occurs at a wave heading near 60° for low 


to moderate seas and at 0° and 180° in severe 
75 
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seas. The horizontal and combined moments, on 
the other hand, exhibit two peaks near UN and 
120 in all sea conditions. For short-crested 
waves, all moment responses show a broad peak 
away from 0° and 180° in low seas. The vertical 
and horizontal moments flatten out in moderate 
seas, and the vertical and combined moment 
assume their biggest values at 0° and 180° in 


Severe Seas. 


The above behavior of the moments with respect to 
wave headings can be explained by the character of 


the RAO. 


The values of the combined moment obtained by the 
two methods described in this paper were not very 
different. Therefore, both methods can be con- 
sidered valid. The first method involves develop- 
ing the RAO for the equivalent combined moment by 
a strip-theory based computational method and per- 
forming the usual power spectral analysis.  Invest- 
igations on the correlation between analytical and 
experimental values have led to the general accep- 
tance of the method first introduced by St. Denis 
and Pierson. However, it is felt that further 
invastigations and refinements are necessary be- 


fore the analytically derived values can be used 
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in actual design. The second method makes use 
of the correlation coefficient, f,,,, to take into 
account the phase difference of the vertical and 
horizontal moments. There is also further need 
to investigate this important statistical 


parameter. 


The comparison study of the predicted and measured 
vertical bending stresses showed good agreement 

in the lower half of the weather spectra en- 
countered in service. In higher seas, the excess- 
ive predicted values are believed to be caused by 
non-linear motions, contrary to the assum»otion of 
the theory. No doubt the addition of bilge keels 
in the analysis of the ship would have reduced the 


CITOIS. 


For the ship's route, both the Roll and Yamanouchi 
corrected wind-to-wave relations seemed to be valid 
representations of the seaway. In view of the ex- 
pected higher-than-actual values by theory, Roll's 
curve is deemed the better approximation of the 
Service route if used in conjuction with the tabu- 


lated data of [24], as was done in this study. 


It was found that there was a significant diff- 


erence between the values of the vertical moment 
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response in long-crested head seas and the 
combined moment response in short-crested seas. 
Again, it is stated that this may not be true 
for other ships, particularly the smaller ones. 
For large tankers with the general character- 
istics of the UNIVERSE IRELAND, the contribution 
of the horizontal moment to the longitudinal 
stresses must be considered in the calculation 


of the primary strength. 


Speed is not a significant factor in the magnitude 


of the wave loads acting on ships. 


It remains to be said that the absolute values of the 
combined moment, or the total longitudinal stress, as pre- 
dicted by the methods in this study could not be compared 
with experimental results. However, the moment trends were 
studied and found to be reasonable and, on the whole, in 
conformity with other researchers' results. Therefore, one 
has confidence on these derived values. The theoretical vert- 
ical bending stress predictions could have been closer to 
measurements but, again, one has firm ideas about the pro- 
bable reasons for the divergence and therefore one is not 


completely dissatisfied. 


The procedure followed in this study is fairly in- 


volved an the cost is not cheap. For further analysis 
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of ships in related studies, it is recommended that the RAO, 
or an approximation, be obtained in the cheapest manner. 
Hopefully, with some anticipation of the expected results, 


more realistic values can be obtained. 
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Appendix A 


SAMPLE COMPUTER OUTPUT DATA 


I. Irregular Wave Results - Long Crested Unidirectional Seas 


Ship Speed = 12 kts. Heading Angle = 22.5 degrees 
Significant Wave Height = 19.2 feet 


Peak Spectral Frequency = 0.5197 


Motions at Origin: 


RMS H(1/3) H(1/10) H(1/1000) 
Heave 0.5974 2,3897 3.0469 4.6003 
Pitch .2456 09822 1.2524 1.8908 
Sway 2237 0.8549 10901 1.6458 
Roll .4386 1.7546 >. 2371 psy 
Yaw 25301 0.5204 076633 70017 
Dynamic Loadings: 
Vertical B.M. 
Sta RMS H (1/3) H(1/10) H(1/1000) 
1 0.0 0.0 Oro D 
2 SII 204.4 260.7 3936 
3 1454.4 58179 7417.9 19599726 
4 8894.8 3557922 45363.5 68490.0 
5 41643.6 166574.5 21 5652.6 3520650670 
6 89928.7 359715 0 458636.6 69245173 
7 144749.5 578998.0 T3822 2.5 TILAS E-O 
8 1° 501.3 799270575 10758987 . 0 15364770 0 
9 25271. 3 100509722 i 0001 1934812.0 
10 29160387 D 1184348.0 1510043.0 22 9669.0 
St 32914 7.6 132865900 Jo7 53.0 25 4457.0 
12 347970.3 139 355170 1774648.0 2673211530 
E 35:9/523. 6 1413294.0 T GST 50.0 2720591: 0 
14 344392.0 1377568.0 155599. 0 26973715. 0 
15 3240996: 7 1284387.0 1637593,0 2472444.0 
16 284297.1 113713820 1449915.0 2189088.0 
17 180498.0 721992.2 920540.1 1396 35. 0 
18 65363.6 261454.4 3338354. 4 50220928 
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Sta. 


T9 
20 
21 
p 


Sta. 


OMANI HAUS WN P= 


Sta ° 


Qn > WW N FP? 


Vertical B.M. 


RMS 


24263.0 
3266.2 
59939: 7 
3083.9 


cont.) 


H(1/3) 


97052733 
1506 .0 
Dyers... 
14735.9 


Horizontal B.M. 


RMS 


0.0 
4568.3 
15469.4 
24915.5 
52746.4 
32.010. 3 
120057. 7 
J 0001 7.7 
201005265 
238483.0 
259256. 3 
209322.5 
269648.6 
260647.3 
242403.2 
214430.1 
132947.2 
48485.1 
22309,9 
6595.0 
456.7 
118.8 


Combined B.M. 


RMS 


0.0 
277.5 
15031.4 
350139.6 
6181.8 


H 1⁄3) 


0.0 

1627277 3.5 
61877.6 
99662.0 
ZEUS GS. 7 
356041.5 
arog STO 
63395130 
834946.2 
953932: 0 
1093025. 5 
1074090.0 
1078594.0 
104258925 
969613.0 
85772025 
53178970 
19394073 
89479.8 
2638052 
1826.93 
475.5 


H(1/3) 


0.0 
16710. 1 
60125.7 

P2052034 
344727.5 


H(1/10) 


12 5741.6 
1692 7.9 
Ze 3. 0 
18788.3 


H(1/10) 


0.0 
21298.8 
78894.0 

1270691 
269006.8 
455952.9 
662784.5 
s 057.5 
1064556.0 
UE D 
J 7107.0 
1369464.0 
1375208.0 
19295010 
J- G256.0 
10935930 
678031.0 
247274.0 
114086.7 
33634.7 
29.2 
606.2 


H(1/10) 


DO 
21305.4 
76660.3 

jn. 0 
01399527. 5 


- 1988574. 


H(1/1000) 


128158 257. 
25100. 
3029. 
28366. 


JONN 


H (1/1000) 


0.0 

E E 
LIS TEA: 
19173829, 
406147. 
68592. 
1000674. 
1316605. 
MEI E. 


19 


— < 


O OOO Odo UI! d On 


20676253. 
2076294.0 
2006984.0 
1866504.0 
1652117.0 
102536953,.7 
57953522 
172248.6 
50781 .8 
3516.7 
91573 


H(1/1000) 


0% 
32167. 
145742. 
23920955 
663600. 


C O F O © 





Combined B.M. 


RMS 


164735.4 
25 106.8 
342644.4 
424373.1 
492267.6 
540433.3 
566890.0 
51411. / 
554399.9 
S1 0609.5 
432582.7 
281527.8 
99230. 1 
37856.3 
( 91.3 
4985.7 

O / 1.7 


(cont. ) 


H (1⁄9) 


658941.7 
1012427.2 
J 05/70 
1697492.0 
1969070.0 
206573370 
220675" 00.0 
228505570 
225599. 0 
2054678.0 
1810331.0 
1l126111.0 

396520.7 

151425.4 

27805.4 
16342.9 
14687.0 


H(1/10) 


840150.7 
1290844.0 
1747486.0 
2164303.0 
2510565.0 
27562 10.0 
2091159.Ü0 
2914210.,0 
2821829.0 
2619712.0 
230810720 
1435792.0 

506073.9 

193067.4 

35451.8 
2083 7.2 
18725.9 


H(1/1000) 


1268462.0 
19489272,0 
2638362.0 
j3267072:0 
3790460.0 
4161336.0 
4365052.0 
4399885.0 
4260409.0 
3955255. 0 
3484887.0 
2167764.0 
764072.3 
291494.0 
53525.4 
31460.1 
282.72.5 
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II. Irregular Wave Results - Short-crested Multi-Directional 
Seas 
Ship Speed = 12 kts. 
Significant Wave Height = 43.3 feet 
Peak Spectral Frequency = 0.3461 
Principal Wind Direction = 22.5 degrees 
Motions at Origin: 
RMS H (15) H(1/10) H(1/1000) 
Heave 6.2544 2520175 31.8974 48.1588 
Pitch ét 4.4734 S@ 1.036 5.656113 
Sway 3,5923 14.3690 7823205 27.6603 
Roll 4.1180 16.4720 2120018 31.7086 
Yaw 06277 225199 322004 4.8334 
Dynamic Loadings: 
Vertical B.M. 
Sta RMS H (1⁄3) H(1/10) H( 171000) 
9 660410.1 2641640.0 336838092. 0 508515820 
10 0714.8 3082859.0 3930646.0 5934504.0 
rl 252958.4 3411833.0 4350088.0 6567779.0 
12 902672.4 43610689 .0 4603629.0 595057 7210 
13 916354.1 3665416.0 4673406.0 205592720 
14 891059.9 356423970 4544406.0 6861161.0 
5 2270132.9 33085310 4219378.0 696822320 
16 728960.0 231532020 2909599965. 0 561299190 
Ly 456989.5 1827358320 2330646.0 35183199 
Horizontal B.M. 
sea RMS H (1⁄3) H(1/10) H(1/1000) 
9 440864.0 1763456.0 2248406.0 3394653.0 
10 geo 02.5 207 610.0 2646403.0 399554970 
11 575644. 3 2302577 <0 2955706. 0 4432461.0 
] 2 606.50.0 2424600.0 3091365.0 46673540 
13 609593.2 2438373.0 3103925. 0 4693867.0 
14 584602.6 2338410.0 2981473.0 4501440.0 
15 333193.7 2140775.0 2529059. O 4120991. 
l6 463863.0 18534520 2365701. 0 3571745. 0 
17 273484.4 109 1917. 0 1394770.0 210583070 





Combined B.M. 


Sta RMS H (1⁄3) Hr 17109 H(1/1000) 
9 958354.0 3833416.0 4887606.0 7579352610 
10 1121474.0 4485896.0 5710551720 8635349.0 
EL 1242414.0 4969656.0 6336 317.0 99650658 7.0 
12 1310208.0 5240832.0 66820610 10088601.0 
13 1320890.0 528356079 Gu 96» 39.0 109870855220 
14 19239550 50937800 6493804.0 9804371. 0 
15 1 / 0845. 0 4683380.0 SSMUS 905506-0 
16 1020343.5 4081374.0 20.752.080 7856645.0 
17 616199.8 2464799.0 3142619.0 4744738.0 





Appendix B 
CALCULATIONS FOR AN AVERAGE VALUE OF er 


For comparison with the value of the correlation 
coefficient of Poy = 0.32 as obtained from figure III.5, 
an average of Poy was calculated from the moment responses 


computed by equation (6) of the first method. 


As the long-term extreme wave moments were not ex- 
trapolated, the mean Qj, value was obtained by averaging 
over wind/wave directions and for representative sea states. 


The results of the calculation are summarized below: 
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V= 16 knots 


Short-crested seas responses 





average m 


mean for low, moderate and high seas 


Pyy = 0.536 





Appendix C 


CALCULATION OF THE MEAN RMS PEAK-TO-PEAK 
VERTICAL BENDING STRESS AND THE 


STANDARD DEVIATION WITH RESPECT TO WIND 


The calculated bending moment was converted into stress 


by the relation: 
Stress = Moment * Section Modulus 


where the section modulus of the ship at the deck is given as 


ola in PBP os 


The computer output gave the RMS of the response. To 
get the peak-to-peak RMS value, the data was multiplied by 


the factor 2/2. 


The mean value of the response to a sea described by 


its significant height is given by: 
= 9 
O CA 
i=l 


where m, = mean RMS of the peak-to-peak response to a sea 


specified by the significant height. 


V. = RMS of the peak-to-peak response to the ith 
> Som spectra 
9 | 
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p(?,)= probability of occurrence of the ith sea spectra. 


p(T. | 1s approximated by the relative frequency of 
occurrence of the different observed wave periods for a 


given wave height in the tabulated data of [24]. 


Average values for all the areas in the ship's route 
were obtained. The calculation is illustrated in table A.l 
for the one case of h = 13.2 feet (wave height code 05 in 


[24]). 


The standard deviation of the vertical bending stress 


with respect to wave measurements is calculated trom the 


relation 
9 
A ~ ~ 
sj = L(v - m.) p(?;) 
1-1 
where s, = standard deviation of the bending stress with 


respect ro waves. 


For direct comparison with experiments based on 


wind measurements, the following relation is used: 


2 
S “ = S 2 + Dë - BH |tan?e [23] 
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The numerical calculation is illustrated below for 


the Roll wind-to-wave relation. 


Beaufort No 


mo wow OQ UN £ t NP 


= 
O 


Beaufort No! h(ft) 


oO N CON C ++ Ç N P 


h(ft) 


m e H H 
00 mn N O 


N) 


O O N O A O O O QV N 


O ON N 
Se oa Oo ` En I 


0 


PQ | ].1712 





2 9 





I 
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123079. 
153520. 
186316. 
229732. 
305014. 
443627. 
541361. 
803314. 
964491. 
1156932. 













.0573 
.0878 
20992 
. 810 

.1184 
. 2064 
. 3912 
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s, (ton-ft) 


- Ui t NY kä - u ON uin 






2 Massey 
E 


9 


. 36x10 


188x107 
.572x10? 
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. 506x10 
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115692 
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-115692 
-115692 
115692 
.115692 
115692 
225692 
11,092 


E 2 
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2 En, - my)" p(T,) 
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in 






.374 
‚414 
.435 
30 
22116 
.668 
. 758 
2 
2:225 
1.403 
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(kpsi) s, (kpsi) 

0.573 .239 
.878 .296 
.0992 315 
.0810 .285 
12184 ‚344 
.2064 ‚454 
39 | a 
-5120 l TO 
. 8637 929 

12.7712 1.032 

s¿(kpsi) s, (kpsi) 

















